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I. Introduction 


The goals of this project are to provide physical insight into the local 
structure of lo's atmosphere, the manner by which gases escape the gravita- 
tional well of the satellite and produce neutral clouds in the circumplanetary 
space, the nature of east-west and possible magnetic longitude asymmetries in 
the plasma torus, and the stability or variability of both the neutral gas 
clouds and the plasma torus over a several year time period. These goals are 
pursued in this AER/JPL collaborative effort by studying and modeling the 
spatial morphology, intensity, and space-time variability of the Io sodium 
cloud as preserved in the JPL Table Mountain Io sodium cloud data set, the 
most complete of the Io sodium cloud data sets currently available. This data 
set documents the 2-D spatial morphology of the D-line emission intensities of 
the cloud on the sky plane as a function of both Io geocetric phase angle and 
the magnetic longitude of Io in the System III coordinate frame over a time 
period from 1976 to 1981. Highly developed and unique models for the Io 
sodium cloud at AER provide the key tool for extracting physical insights from 
this data set. 

The complete JPL Io sodium cloud data set is composed of 1974-1979 
spectral data, 1981 spectral data, 1976-1979 Region B/C image data, 1981 
Region B/C image data, and 1981 Region A image data. The complete data set 
was reviewed and assessed in the first project year (see the 1985 Annual 
Report). From this review, both the quality and quantity of data acquired in 
1981 were in general much superior to the earlier observations. Because of 
this fact, the primary emphasis of modeling analysis in this project has been 
directed toward the 1981 data set and, in particular, to the 1981 Region B/C 
images. Efforts in the second year have therefore focused primarily on the 
1981 Region B/C image data. 

Progress in the second project year can be divided into three main 
areas: (1) data quality review for the 1981 Region B/C images, (2) selection 

and data processing of a subset of 1981 Region B/C images, and (3) physical 
enhancement of the Io sodium cloud model in preparation for a third year 
analysis of the data. A major step in the second area was the successful 
development and application of software at the Multimission Image Processing 
Laboratory (MIPL) of JPL to remove distortions present in the brightness 



distribution of each image. Because there were significant and unavoidable 
delays encountered in developing this image software in the past two years, 
the project was extended for a third year to provide the necessary time to 
analyze the processed image data. The following more detailed discussion of 
second year progress is divided into the three main areas noted above. 



II. Data Quality Review in the 1981 Region B/C Images 

Observational parameters for the 263 Region B/C images of the 1981 JPL 
data acquired over 14 nights are summarized in numerical form in Table 1 and 
in graphical form in Figure 1. Prior to this project most of these images had 
been processed only preliminarily to remove background signals. Images on May 
4, 5, 12, 13 and June 4 had undergone further processing to remove the 
instumental response function and to normalize the image intensity. Because 
of this, the 110 images (42% of the data set) acquired on May 4, 5, 12 and 13 
were chosen and thoroughly reviewed during the first project year. As a 
result of this review, further refinements in the image processing of these 
data frames were indentified. The quality of the remaining 153 images is 
reviewed here as a prerequisite for selection of images in our model-data 
studies. 

The results of this data quality review are summarized in Table 2 and 
were compiled by manual inspection of photographical products for each image 
produced prior to this project by the Image Processing Laboratory of JPL. 

Most images are composed by adding three data frames, where each data frame 
generally represents an exposure time of approximate ten minutes. Some images 
are, however, composed of only one or two data frames. The quality of the 
images is rated in the right hand column of Table 2 and qualifying comments 
are also indicated. As can be seen, almost all of the images are of 
sufficient quality to be valuable in our analysis of this data set. Those 
images selected for more detailed studies will, as noted earlier, undergo 
further image processing. 
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1981 Region B/C Images: Observing Chronology 
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System III ( 1965) 


System III Magnetic Longitude of lo 


1981 10 SODIUM CLOUD DATA 
FROM TABLE MOUNTAIN OBSERVATORY 



lo Geocentric Phase Angle 


Figure 1. Observing Parameters for 1981 lo Sodium Cloud Images. 

The angular coverage for the lo geocentric phase angle and the System III 
magnetic longitude of lo over which Region B/C images were recorded In the 
JPL Table Mountain Data Set is indicated for all 14 nights of observations. 




Table 2 

1981 Region B/C Inages: Data Quality Review 


Date (UT) 

UT Start 

UT End 

Frames 

Quality 

March 25 

- 

- 

— 

_ 1 


5:11 

8:15 

29 

usable^ 


_ 



_ 1 


*only raw data available; 16 images with most 

exposures 1 1/2 minutes. 

since images were so 

close to Jupiter; 

must 

add images for useful 

results; images were 

heavily masked in 

both 

east and west directions. 

April 5 5:18 

5:22 

18 

O 

marginal 

5:25 

5:32 

19 

2 

good 

5:36 

5:46 

20 

2 

good 

5:48 

6:03 

21 

2 

good 

6:05 

6:20 

22 

good (some saturation) 

6:24 

6:44 

23 

good (some saturation) 

6:47 

7:02 

24 

2 

good 

7:04 

7:14 

25 

2 

very good 

7:17 

7:24 

26 

very good^ 

7:26 

7:36 

27 

very good^ 

7:40 

7:50 

28 

2 

very good 

O 

Part of forward cloud 

truncated by mask. 


April 6 (one) 4:33 

4:40 

11 

good 

4:42 

4:52 

12 

good"* 

4:54 

5:09 

13 

3 

good 

5:11 

5:28 

14 

3 

good 

5:30 

5:50 

15 

good"* 

5:52 

6:09 

16 

3 

good 

(two) 6:38 

6:55 

22 

4 

good J * 

6:59 

7:09 

23 

good^ 

7:11 

7:23 

24 

good^ 

7:25 

7:42 

25 

good^’^ 

7:44 

8:01 

26 

good^ 


J only raw data available. 

^probably part of forward cloud truncated 



Table 2 
(continued) 



Date (UT) 

UT 

Start 

UT End 


Frames 

Quality 


April 28 


4:21 

5:01 


25-27 

good 

• 



4:37 

5:13 


26-28 

very good 




4:49 

5:24 


27-29 

very good 




5:03 

5:34 


28-30 

very good 




5:15 

5:45 


29-31 

very good 




5:25 

5:56 


30-32 

very good 




5:36 

6:07 


31-33 

very good 




5:47 

6:18 


32-34 

very good 




5:58 

6:30 


33-35 

very good 




6:09 

6:40 


34-36 

excellent 




6:31 

6:52 


36, 38 5 

very good 




6:43 

7:03 


38, 39 6 

very good 




6:54 

7:16 


39, 41 7 

very good 




7:07 

7:27 


41, 42 6 

excellent 




7:07 

7:37 


41-43 

excellent 




7:18 

7:46. 

5 

42-44 

very good 




7:28 

7:57. 

5 

43-45 

very good 




7:39 

8:07 


44-46 

very good 




7:49 

8:18 


45-47 

very good 




7:59 

8:28 


46-48 

good 




8:10 

8:38 


47-49 

good 




8:20 

8:48 


48-50 

good 




8:30 

8:58 


49-51 

good 


"’two 

frames 

only with 

9 minute 

gap 

in the middle. 



6 two 

frames 

only. 






7 two 

frames 

with 1 minute gap in center. 



April 29 


3:18 

3:43 


1-3 

usable 




3:26 

3:57 


2-4 

good 




3:35 

4:07 


3-5 

good 




3:48 

4:19 


4-6 

good 




3:58 

4:32 


5-7 

very good 




4:10 

4:42 


6-8 

very good 




4:22 

4:53 


7-9 

very good 




4:32 

5:05 


8-10 

very good 




4:44 

5:16 


9-11 

very good 1 




4:55 

5:32 


10-12 

very good 




5:06 

5:43 


11-13 

very good 




5:21 

5:55 


12-14 

very good 

• 



5:33 

6:08 


13-15 

very good 


7 



Table 2 
(continued) 


Date (UT) UT Start 

UT End 

Frames 

Qual 

April 29 (cont.) 5:45 

6:21 

14-16 

good 

5:57 

6:32 

15-17 

good 

6:10 

6:42 

16-18 

good 

6:23 

6:52 

17-19 

good 

6:33 

7:02 

18-20 

good 

6:43 

7:14 

19-21 

good 

6:53 

7:25 

20-22 

good 


“background level increasing. 

^little portion of forward cloud truncated by mask. 


May 6 


3:07 

3:25 

5-7 

usable (twilight) 

3:12 

3:36 

6-8 

good 

3:18 

3:46 

7-9 

good 

3:27 

3:57 

8-10 

very good 

3:37 

4:08 

9-11 

excellent 

3:47 

4:20 

10-12 

excellent 

3:58 

4:32(?) 

11-13 

excellent 

4:09 

4:44 

12, 13 10 

good 

4:22 

5:01 

13, 15 10 

good 

4:50 

5:12 

15, 16 10 

good 

4:50 

5:24 

15-17 

excellent 

5:03 

5:35 

16-18 

excellent 

5:14 

5:47 

17-19 

excellent 

5:25 

5:59 

18-20 

excellent 

5:36 

6:10 

19-21 

excellent 

5:48 

6:20 

20-22 

excellent 

6:00 

6:31 

21-23 

excellent 

6:11 

6:43 

22-24 

excellent 

6:22 

6:54 

23-25 

excellent 

6:33 

7:04 

24-26 

excellent 

6:44 

7:14 

25-27 

very good 

6:55 

7:28 

26-28 

very good 

7:05 

7:40 

27-29 

very good 

7:18 

7:52 

28-30 

very good 

7:30 

8:05 

29-31 

good 

7:41 

8:16 

30-32 

good 

7:54 

8:26 

31-33 

usable^ * 

- 

- 

32-34 

no good 


^no signal in frame 14. 
Jupiter close to horizon. 


8 





Table 2 



(continued) 


Date (UT) UT Start 

UT End 

Frames 

Quality 




12 

June 4 “ 




4:57 

5:35 

12-14 

usable^ 

5:09 

5:47 

13-15 

13 

usable 

5:24 

6:00 

14-16 

usable^ 

^perhaps more images. 

but must 

examine date 

on tapes. 

1 3 

forward cloud may be 

partially truncated. 


June 5 (one) 3:38 

4:05 

5-7 

marginal^ 

3:47 

4:18 

6-8 

good^ 

(two) 5:27 

6:09 

11-13 

good 1 "* 

5:40 

6:25 

12-14 

good 

5:56 

6:39 

13-15 

good 

6:12 

6:54 

14-16 

good 

6:26 

7:07 

15-17 

good 

Idiots of background noise. 



15 some background noise 

1 • 



June 6 

- 

- 

_ 16 

4:08 

4:41 

9-11 

very good 

4:19 

4:53 

10-12 

very good 

4:30 

5:05 

11-13 

very good 

4:42 

5:23 

12-14 

very good 

4:54 

5:42 

13-15 

very good 

5:08 

5:55 

14-16 

very good 

5:24 

6:07 

15-17 

very good 

^possibly more images. 




June 14 3:38 

4:04 

5-7 

usable^ » 

3:47 

4:15 

6-8 

good^ 

3:56 

4:28 

7-9 

good 

4:05 

4:40 

8-10 

good 

4:17 

4:52 

9-11 

good 

4:29 

5:05 

10-12 

good 

4:41 

5:19 

11-13 

good 

4:53 

5:32 

12-14 

good 


^signal to noise not very good 

^portion of forward cloud truncated by mask 



III. Data Processing Activities 

1 . Data Correlation Studies 

Five data correlation studies are defined in Table 3 for the set of 1981 
Region B/C images. Each study provides a different classification for these 
images based upon their Io geocentric phase angle, 4>, and their Io System III 
magnetic longitude angle, i|>. The five classifications are useful in identi- 
fying and then studying the effects of radiation pressure, the atom source 
characteristics, and the plasma torus sink characteristics on the morphology 
of the sodium cloud. Each correlation study lists those images that are 
appropriate for the defined classification, regardless of the image data- 
quality summarized in Table 2. Incorporating the results of Table 2 will then 
aid in selection and prioritization of images to undergo further image 
processing and then model-data analysis. 

The results of the five data correlation studies of Table 3 are 
summarized in Tables 4-8, respectively. From Table 4, consecutive image 
studies are possible on all 14 nights, but short exposure times on March 25 
and questionable data quality on June 4 might eliminate these two evenings 
(see Table 2). From Table 5, cloud stability studies are available for three 
image sequence pairs (also see Figure 1) which all have good to excellent data 
quality. If comparisons of these image sequence pairs indicate that the 
sodium cloud and its interaction with the plasma torus are stable over a few 
month time interval, then the procedure of comparing images of different dates 
in the correlation studies III, IV and V will represent a valid multi-faceted 
look at a common (although periodically time dependent) cloud-torus state. 

From Table 6, east-west image studies of the cloud morphology may be 
undertaken and used to explore the nature of east-west difference in the 
plasma torus also identified by optical torus emissions (Morgan and Fertel, 
1984; Morgan, 1985) and by UV emission from the Io plasma torus (Sandel and 
Broadfoot, 1982; Shemansky and Sandel, 1982). From Table 7, System III 
variations in the electron properties of the plasma torus may be investigated 
by comparison and analysis of images. From Table 8, the effects of viewing 
angle on the cloud morphology may be investigated. The effects of system III 
variation can also be investigated by comparing images for somewhat similar 
phase angles, and excellent cases exist for <p in the range from about 160° to 
310° (see Table 8 and Figure 1). 



Table 3 


1981 Region B/C Inages: Five Correlation Studiest 



I. 

Consecutive Image Study: 

<|>(t) , y(t) 

• 

II. 

Cloud Stability Study: 4^) * ( <J> 2 » 4^) 

image sequences separated by 32 days 


III. 

East/West Image Study: 

♦ w = <fr E + 180° 

• 


^E - *W 

V w - * E + 180° 

only one image-sequence pair 
no image pairs 



^E * 

many comparison images 


IV. System III Variability Study: <j>^ = <J> 2 » 4^ * V ^ 

V. Io Phase Angle Variability Study: $ 1 * < l > 2» = ^2 


t<j> = Io Geocentric Phase Angle 
V = Io System III Magnetic Longitude Angle 




11 



Table 4 

1981 Region B/C Images: Correlation Study I 


Consecutive Image Study : <t>(t), t ) 

Image data suitable on all 14 nights 


1 . 

March 25 



2. 

April 5 



3. 

April 6 



4. 

April 28 



5. 

April 29 



6. 

May 4 

(June 

5) 

7. 

May 5 

(June 

6) 

8. 

May 6 



9. 

May 12 



10. 

May 13 

(June 

14) 

11. 

June 4 





Table 5 

1981 Region B/C Inages : Correlation Study 11 


Cloud Stability Study : i|)j) “ ^2^ 


Three Image Sequence Pairs Suitable 


(May 5, June 6) 
(May 4, June 5) 
(May 13, June 14) 


Io east of Jupiter 
Io west of Jupiter 


Io west of Jupiter 



Table 6 

1981 Region B/C Images : Correlation Study III 

East/West Image Study : 4>w = ‘t’E + 


East Images 
<)>£ 

March 25 
May 12 
June 4 
May 5 
June 6 
April 28 


West Images 

^E 
May 4 

June 5 (one) 
June 5 (two) 
May 13 
May 6 

April 6 (one) 
April 6 (two) 
June 14 
April 29 


♦w 

May 4 

May 4, May 6, May 13, June 5 (one, two) 

May 4 

April 29, May 4, May 6, May 13, June 5 (one, two) 
April 6 (one), April 29, May 4, May 13, June 14 
April 6 (one)*, April 29, May 4, May 6, May 13, 
June 5 (two), June 14 


March 25, April 5, April 28, May 5, May 12, June 4 
May 5, May 12 

April 5, April 28, May 5, May 12 

April 5, April 28, May 5, May 12, June 6 

April 5, April 28, May 5, May 12, June 6 

April 5, April 28*, May 5, June 6 
no candidates 

April 5, April 28, May 5, June 6 

April 5, April 28, May 5, June 6 


*iJj e « in addition 



Table 7 

1981 Region B/C Inages: Correlation Study IV 


System III Variability Studyt : 4*1 = <}>2* ^1 * ^2 


Io East of Jupiter 


Io West of Jupiter 


4 > 

Date 


4 > 

Date . . ¥ 

50 

June 4 

26 

250 

May 

13. .102 


May 12 

167 


May 

4. .323 

60 

June 4 

58 

255 

May 

13. .118 


May 12 

198 


May 

4. .338 





June 

5. .338 

65 

June 4 

72 





May 5 

172 

260 

May 

13. .133 


May 12 

212 


May 

4. .353 

70 

May 5 

190 

270 

June 

5. ..25 


May 12 

231 


May 

4. ..25 





May 

6. .128 

80 

May 5 

222 


May 

13. .167 


May 12 

263 

280 

May 

6. .157 

90 

April 5 

5 


May 

13. .197 


April 28 

219 





May 5 

256 

290 

April 

29. .153 


May 12 

294 


May 

6. .190 





May 

13. .232 

100 

April 5 

37 


June 

14. .232 


April 28 

251 


April 

6 (one) 


May 5 

288 





June 6 

288 

300 

April 

29. .185 





May 

6. .223 

110 

April 5 

70 





April 28 

283 

305 

April 

29. .202 


May 5 

322 


May 6.. 238 


June 6 

322 


April 

6(two) 




310 

April 

6. . . .4 





April 

29. .219 


299 


350 


tangles <|> and given in degrees 
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Table 8 


1981 Region B/C Images : Correlation Study V 


Io Phase Angle Variability Studyt : <J>i * <f> 2 » * $2 


4* 

Date 

♦ 

1 f 

Date .. 4) 







0 

May 4 

262 

140 

May 13 

262 


April 16 (two) 

309 


May 6 

274 

10 

May 4 

265 


April 29 

286 


April 6(two) 

312 

150 

May 12 

45 

20 

June 4 

48 


May 13 

265 


April 5 

94 f 


May 6 

277 


May 4 

268 f 


April 29 

289 


June 5(two) 

268? 

160 

May 12 

48 

30 

June 4 

51 


May 5 

60 


April 5 

97l 


May 13 

268 


May 4 

27lf 


May 6 

280 

40 




April 29 

292 

June 4 

54 




April 5 

100 

170 

May 12 

51 

50 




May 5 

63 

June 4 

57 


271 


May 13 


April 5 

103 






May 6 

283 

60 

June 4 

60 


April 29 

295 


April 5 

106 

180 

May 12 

54 

70 

June 4 

63 


May 5 

66 


April 5 

109 


May 13 

274 

110 

May 13 

253 


May 6 

286 


May 6 

266 


April 29 

298 

120 

May 13 

256 

190 

May 12 

57 


May 6 

268 


May 5 

69 


April 29 

280 


May 13 

277 

130 




May 6 

289 

May 13 

259 






April 29 

301 


May 6 

271 





April 29 

283 






Table 8 
(continued) 




200 


210 


220 


230 


240 


Date 


May 12 

60 

May 5 

72 

May 13 

280 

May 6 

293 

April 29 

304 

May 12 

63 

May 5 

75 

April 28 

87 

May 4 

216 

May 13 

283 

June 14 

283 

May 6 

296 

April 29 

307 

May 12 

66 

May 5 

78 

April 28 

90 

May 4 

219 

May 13 

286 

June 14 

286 

May 6 

299 

April 29 

310 

May 12 

69 

May 5 

81 

April 28 

93 

May 4 

222 

May 13 

289 

June 14 

289 

May 6 

302 

May 12 

72 

May 5 

84 

April 28 

96 

May 4 

225 

May 13 

292 

June 14 

292 

May 6 

305 


¥ 

Date .. <t> 


250 

May 12 

75 


May 5 

87 


April 28 

99 


May 4 

228 


June 14 

295 

260 

May 12 

78 


May 5 

90 


April 28 

102 


May 4 

231 

270 

May 12 

81 


May 5 

94 


April 28 

105 


May 4 

234 

280 

March 25 

26 


May 12 

84 


June 6 

96 


May 5 

97 


April 28 

108 


May 4 

237 

290 

March 25 

29 


May 12 

87 


June 6 

99 


May 5 

100 


April 28 

111 


May 4 

240 

300 

March 25 

32 


May 12 

90 


June 6 

102 


May 5 

103 


April 28 

114? 


May 4 

243 


April 6(one) 

29ol 



Table 8 
(continued) 



l 

Date 

♦ 

1 * 

Date . . <}> 



310 

March 25 

36 

330 

June 6 

111 

A 


June 6 

105 


May 4 

252 

w 


May 5 

106 


June 5(one) 

252 



April 28 

117? 

350 

May 4 

258 



May 4 

246 


April 6(two) 

306 

• 


April 6(one) 

294l 





320 

June 6 

108 






May 5 

109 






April 28 

120$ 




• 


May 4 

249 






April 6(one) 

297$ 





tAngles <j> and ij> given In degrees 

flo phase angles differ by approximately 180° 
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2. Development of Data Processing Software 

Data processing priorities defined for the JPL Table Mountain Io Sodium 
cloud data are summarized in Table 9. The first and second priorities require 
development of special software at MIPL. 

The development of software for removal of distortion in the D2 and D^ 
image data (item 1 of Table 9) was successfully completed in the third 
quarter. This completion marked a major milestone for the project and was by 
far the most difficult task of Table 9* This development had proceeded much 
slower than initially expected, partly because of its complexity, but primary 
because of delays at MIPL beyond our control. These delays were caused by two 
main factors: (1) major changes in the MIPL computer hardware (and corres- 

ponding conversion of all of the MIPL image processing software) that occurred 
in the first year and a half of this project, and (2) difficulties in schedul- 
ing time of sufficiently senior MIPL analysts for actual software develop- 
ment. The new image software is quite complex, and its application to image 
data requires interactive analyst intervention for each image processed. 

The development of software to determine the absolute brightness cali- 
bration of an image (item 2 of Table 9) was initiated near the end of the 
third quarter and completed in the beginning of the fourth quarter. The 
procedure is straight forward and has been thus far restricted to the 1-D slit 
data which is sufficient for our calibration process. 

3. Application of Data Processing Software 

The third data processing priority for the JPL Table Mountain Io sodium 
cloud data in Table 9 is the application of the image and calibration software 
to images in preparation for modeling analysis. This application began in the 
fourth quarter at MIPL and was restricted to 1981 Region B/C images data. At 
the end of the second project year, twenty-two of these image have been 
completely processed and several additional images are in progress. The 
twenty-two completed images are identified in Table 10 in terms of their 
observation date, integration start and end times, and their SIP tape number 
plus the ID number of each frame that is co-added to produce the image. The 
start and end values of the geocentric Io phase angle and system III longitude 
angle of Io are also given in Table 10. The twenty-two images were selected 
to yield a useful subset of images appropriate for the data correlation 
studies discussed above. The absolute brightness calibration for each of these 
images and others currently being processed are expected to be available early 
in the third project year. 


Table 9 

Data Processing Priorities 


Development of Image Processing Software to Remove Distortion in the 
Brightness Distribution of Images 

a. Implement improved background subtraction techniques for images. 

b. Improve techniques for removal of image distortion near Io 
produced by continuum light scattered by Io. 


Development of Software to Determine the Absolute Brightness Calibration 
of Images 

a. Analyze 1981 1-D slit data on Io's disk and Region A image data to 
establish an absolute brightness calibration for the 1981 data. 

b. Analyze 1976-79 1-D slit data on Io's disk to establish a 
brightness calibration for the 1976-79 data. 


Application of Software to Prepare Images for Modeling Analysis 

a. Using improved techniques in 1 above, remove brightness morphology 
distortions in a selected subset of images. 

b. Using the information in 2 above, absolutely calibrate this 
selected subset of images. 
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IV. Modeling Activities 


1 . Improvements for the To Sodium Cloud Model 

The lifetime description of sodium in the planetary magnetosphere has 
been improved in the second year to include the offset dipole nature of the 
magnetic field as well as the effects of the east-west electric field (Barbosa 
and Kivelson 1983; Ip and Goertz 1983). Inclusion of these factors in the 
lifetime description is important to model properly the morphology of the 
sodium cloud as discussed below. 

The offset magnetic dipole causes the plasma torus center to be offset 
from the center of Jupiter. This causes the lifetime of sodium at Io's 
position to have a slight (but true) system III modulation in addition to the 
system III oscillation caused by the ~7° tilt between the rotating centrifugal 
plasma plane and the satellite orbit plane. The lifetime variations due to 
this oscillation were discussed earlier (see third progress report of the 
first year). The lifetime modulations due to the dipole offset are relatively 
small compared to the effects produced by the oscillation. These modulations 
will modify the amount of sodium present in the cloud but not its basic 
morphology. 

The east-west electric field has three major effects on the plasma 
properties: (I) a geometric shift of the plasma torus center toward the east, 

(2) a decrease in the electron temperature of a plasma volume element as it 
moves from west to east of Jupiter, and (3) a decrease in the electron number 
density of a plasma volume element as it moves from west to east of Jupiter. 
These three effects provide a significantly longer lifetime for sodium atoms 
near Io when the satellite is east of Jupiter than when the satellite is west 
of Jupiter. This east— west asymmetry in the sodium lifetime (because of 
electron impact ionization) provides a natural explanation for the east-west 
intensity asymmetry of the sodium cloud that was initially discovered by 
Bergstralh (1975, 1977) and noted more recently (in a more spatially extended 
region about Io) in the JPL Table Mountain sodium cloud image data (Goldberg, 
Garneau and LaVoie 1984). 


2. Preparation for Modeling Analysis In the Third Year 

In preparation for third year analysis of image data, model calculations 
have been performed at AER in the third and fourth quarters of this project 
year to examine the impact of the new lifetime description (discussed above) 
on the intensity and morphology of the sodium cloud. To describe the model 
and the nature of the space/time dependent lifetime on the cloud, two 
companion papers have been undertaken. The first paper entitled "A General 
Model for Io's Neutral Gas Clouds: Mathematical Description" is a general 

model-documentation paper. This paper has been completed and is included in 
the Appendix. The second paper entitled "A General Model for Io's Neutral Gas 
Cloud: Application to the Sodium Cloud" explicitly applies the general cloud 

model to the sodium cloud. The model calculations simulate the basic cloud 
morphology and the character of the east-west orbital and east-west intensity 
asymmetries for the sodium cloud that have been previously reported in the 
literature. The second paper is approximately 80% completed and will be 
finished in the third project year at which time these two companion papers 
will be jointly submitted for publication. These publications will extend the 
more recent AER theoretical contributions to the understanding of the Io 
sodium cloud (Smyth 1979, 1983; Pilcher et al., 1984; Smyth and Combi 1983, 
1984, 1985) and provide a solid basis for the third year analysis. 
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APPENDIX 


A GENERAL MODEL FOR 10' S NEUTRAL GAS CLOUD 
MATHEMATICAL DESCRIPTION 



A General Model for Io’s Neutral Gas Clouds: 
Mathematical Description 


by 


W • H. Smyth 
M. R. Combi 


Atmospheric and Environmental Research, Inc. 
840 Memorial Drive, Cambridge, Massachusetts 


[PRELIMINARY DRAFT] 


ABSTRACT 


A general mathematical formalism for calculating the physical properties 
of any of Io's neutral gas clouds (Na, K, 0, S, S0 2 ) i S presented. The 
dynamical effects of both the gravitational fields of Io and Jupiter and of 
solar radiation pressure are included in addition to the many complex space 
and time dependent interactions that occur between the neutral clouds and the 
plasma torus. Earlier models have not included these complex space and time 
dependent interactions. The importance of this new model in studying both the 
plasma conditions prevalent in the inner planetary magnetosphere and the 
nature of Io's local atmosphere is discussed. A numerical method for 
evaluating the physical properties of the neutral clouds has been developed 
and is described. The mathematical formalism and numerical method are applied 
to a study of the Io sodium cloud in a companion paper in this issue. 



« 

1 . Introduction 

During the past two decades, Io, the innermost of the four Galilean 
satellites of Jupiter, and its profound influence on the planetary magneto- 4 

sphere have been one of the most exciting and intensely studied subjects in 
the solar system. The first evidence that Jupiter and Io constituted a com- 
plex interactive system was introduced by Bigg (1964), who discovered that 

some of the decametric radio emission of Jupiter was directly influenced by 4 

the satellite’s orbital and magnetic location about the planet. Since then, a 

large information base has been accumulated from ground-based observations, 

rocket measurements, earth— orbiting satellites, and four spacecraft encounters 

with Jupiter: Pioneer 10 (4 December 1973), Pioneer 11 (2 December 1974), 4 

Voyager 1 (5 March 1979), and Voyager 2 (9 July 1979). 

t 

From ensuing studies of this information base, a picture for the Io- 
Jupiter system is emerging and slowly coming into focus. Io has a local atmo- 
sphere, although it is not well characterized or structurally understood. Io 4 

has extended neutral atmospheres or clouds that exist large distances beyond 
the gravitational grasp of the satellite and fill partial or somewhat incom- 
plete toroidal-shaped volumes around Jupiter. These clouds occupy an overlap- 
ping volume with the Io plasma torus - the most dense part of the planetary 4 

magnetosphere - and interact with it through collisional processes including 
ionization and charge exchange reactions. These interactions not only shape 
the spatial morphology of the clouds but supply heavy ions to the torus and 
affect the composition, structure, and energy budget of the plasma torus and 4 

its extension into the larger magnetosphere. The primary heavy ions that com- 
pose the corotating (or nearly corotating) plasma torus are known to be oxygen 
(0 + , 0 ++ ) and sulfur (S + , S ++ , S +++ ) , although the relative abundances of 

these ions throughout the torus volume (even for a two-dimensional longitudin- 4 

ally-averaged picture) are still not well established. Io and its local atmo- 
sphere also interact with the corotating plasma torus and planetary magnetic 
field in ways that are yet to be fully understood. Some of these interactions 
almost certainly provide ions and energy directly to the plasma torus as well f 

as drive escape mechanisms for gases in the local atmosphere. 

For purposes of this paper, a brief review of Io's local and extended 
atmospheres and their relevant interactions with the plasma torus will be 

given in section 2. This will provide necessary information for understanding ( 

the evolution and motivation for developing the new general neutral gas cloud 
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model presented in sections 3 and 4. A more complete discussion of lo may be 
found in several excellent review articles (Nash et al. 1986; Brown, Pilcher 
and Strobel 1983; Fanale et al. 1982, 1977; Kumar and Hunten 1982; Pilcher and 
Strobel 1982; Brown and Yung 1976). The paper concludes with a brief summary 
in Section 5. The hew model is applied to the sodium cloud in a companion 
paper (Smyth and Combi 1986) in this issue. In a series of future papers, the 
new model will also be applied to other neutral clouds of lo (K, 0, S, SO 2 ). 



2. Historical Perspective 

Our understanding of Io's local atmosphere, although far from complete, 
has grown steadily since 1964, when Binder and Cruikshank (1964) reported an 
anomalous brightening of the satellite's surface as it emerged from eclipse. 
This post-eclipse brightening suggested a dynamic and volatile atmosphere 
for Io, for which in 1971 an upper limit surface pressure of approximately 
10 - ^ bars was deduced by Smith and Smith (1972) from occultation of the star 
S Scorpii by the satellite. In the summer of 1972, the discovery by Brown 
(1974) of sodium optical emission from Io's vicinity and the subsequent obser- 
vations in 1973 by Trafton, Parkinson and Macy (1974) of an extended sodium 
cloud well beyond the gravitational reach of Io confirmed the presence of a 
local atmosphere. This discovery also sparked a large number of earth— based 
observations of the sodium cloud that have continued over the intervening 
fourteen years. An ionsophere for Io discovered from analyses of Pioneer 10 
S-band radio occultation data (Kilore et al. 1974) provided additional 
information to characterize the satellite's local atmosphere and its possible 
interactions with the Jovian magnetosphere (Cloutier et al. 1978). A quantum 
leap forward in understanding Io's local atmosphere was obtained by the 
Voyager 1 spacecraft's discovery of active volcanic plumes (Morabito et al. 
1979) and detection of SO 2 gas (~ 10 - ^ bars surface pressure) over one of many 
hot spots on Io (Pearl et al. 1979). Although SO 2 now appears to be the major 
constituent of the atmosphere (Kumar 1979), the vertical structure, global 
extent and surface density variations of the atmosphere still remain uncertain 
even after much theoretical analysis (Kumar 1980, 1982, 1984, 1985; Summers 
1984; Matson and Nash 1983). Insightful information for the density of sodium 
in Io's local atmosphere is to be expected in the near future, however, upon 
analysis of novel earth-based observations recently presented by Schneider 
et al. (1985) of absorption signatures at the sodium D-line wavelengths of 
reflected sunlight from Europa as it was eclipsed by Io. 

The discovery of the Io sodium cloud marked the first of four extended 
atmospheres that have been detected to date for the satellite as summarized in 
Table 1. The three more recently detected neutral clouds for Io are potassium 
in 1975 and atomic oxygen in 1980, both from ground-based observations of 
optical lines, and atomic sulfur in 1981 from rocket-borne instruments 
measuring ultraviolet emissions. A fifth cloud of neutral SO 2 * s expected, 
because of the presence of a small amount of S 02 + ions measured in the plasma 



torus (Bagenal and Sullivan 1981; Bagenal 1985), but has not been detected 
directly. For sodium, more than 1000 observations of its spatial brightness 
have been made in the last 14 years, with most of the earlier measurements 
providing only slit-averaged intensity values and most of the more recent 
measurements providing increasingly higher-quality two-dimensional image 
data. For the dimmer potassium cloud (~ 1/20 the brightness of the sodium 
cloud), about 200 slit-averaged observations have been made almost entirely by 
Trafton (1975, 1981), who presented a survey of this gas cloud that exhibited 
spatial and temporal variations similar to those documented for the Io sodium 
cloud. For the much dimmer atomic oxygen and atomic sulfur clouds, only a 
very few slit-averaged detection observations have been acquired. 

The presence of neutral clouds for Io provides a largely unexplored but 
rich potential for studying both the plasma conditions prevalent in the inner 
planetary magnetosphere and the nature of Io's local atmosphere. Sodium is by 
far the best neutral cloud to use as a probe for these purposes, because it is 
the brightest and most extensively observed and is (together with potassium) 
sensitive to spatial and temporal changes in the plasma torus due to its small 
electron impact ionization lifetime (~ 1-2 hours) near Io's orbit. Early ob- 
servations acquired in 1973 and 1974 (i.e., pre- and post-Pioneer 10) by Macy 
and Trafton (1975 a,b) and Mttnch and Bergstralh (1977), for example, first 
suggested the presence of a magnetospheric sink for sodium near Io's orbit 
that was both radially asymmetric and oscillatory. The radial asymmetry was 
suggested by the presence of a predominatly forward sodium cloud, and the 
oscillatory nature was suggested by time dependent changes in the sodium 
brightness north and south of the satellite that were correlated with the 
satellite's magnetic latitude (Trafton and Macy 1975; Trafton 1980). As will 
be shown in the companion paper in this issue, the predominantly forward cloud 
is indeed naturally produced by the radial asymmetry of the plasma torus for 
sodium ejected from Io at relatively slow speeds (~ 3 km sec - ^) as first sug- 
gested from these early observations. These slow speed atoms may be energized 
by magnetospheric ions both directly through scattering of neutrals from the 
local atmosphere (or surface, if regions of low density exist) as discussed in 
a number of papers (Sieveka and Johnson 1984; Cheng 1984; Lanzerotti and Brown 
1983; Brown, Pilcher and Strobel 1983; Ip 1982; Matson, Johnson and Fanale 
1974), and indirectly through multicollisional processes resulting in atmo- 
spheric sputtering (Summers, Yung and Haff 1983; Haff, Watson and Yung 1981; 



McGrath and Johnson 1986) or blowoff (Hunten 1985). More recent observations 
of the sodium cloud have shown , for example, (1) well-defined spatial zones 
that become temporally deficient of atoms, due perhaps to locally (likely 
System III) enhanced ionization (Goldberg, Gameau, and LaVoie 1984), (2) 
directional features produced by high-velocity (~ 20 km sec - *) sodium driven 
elastically by the magnetospheric wind from the local atmosphere at near right 
angles to lo's orbit and preferentially from the equatorial regions (Pilcher 
et al. 1984), and (3) very fast (£ 60 km sec -1 ) narrowly-collimated jets of 
sodium driven directly ahead of the satellite by charge exchange reactions 
with the plasma torus (Trauger 1984, 1985). Although the oxygen and sulfur 
clouds are only bright enough to be barely detected, physical insight acquired 
for the bright sodium cloud may be adopted where appropriate in 01, SI and SO 2 
cloud models. These atomic and molecular cloud models so constructued then 
become enhanced tools themselves to determine ion loading rates, plasma mass 
loading rates, and ion energy input rates to the plasma torus as well as to 
specify the character and source rate of neutral gases escaping from lo's 
local atmosphere (Smyth and Shemansky 1983). 

Models for the neutral clouds of lo are essential in probing and 
recovering physical information from spatial brightness observations of lo's 
extended atmospheres and relating this information to conditions in lo's local 
atmosphere and the planetary magnetosphere. The development of cloud models 
for Io was historically initiated for the sodium cloud in the early 1970' s and 
has evolved during the last decade with progressively more complete descrip- 
tions for the source, lifetime, and orbital dynamics of sodium atoms in the 
Jupiter system. Only in the last few years have models for the Io oxygen and 
sulfur clouds (Smyth and Shemansky 1983) and the Io potassium cloud (Smyth and 
Combi 1984) been developed. A general model, appropriate to any of these neu- 
tral gas clouds of Io, has now been formulated and is presented in Sections 3 
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Mathematical Description for the lo Neutral Gas Clouds 

Mathematical Approach 

The basic mathematical approach is based upon the conservation in phase 

space of the one-particle distribution function f(v,x,t) for a gas species. 

In brief, if at an initial time t = t, the distribution function can be 

b 

specified, then the distribution function at a later time may, under certain 
circumstances of interest here, be simply related to its initial value by 
properly mapping its time evolution in phase space and including processes 
that modify its initial value during this time evolution. This conservation 
in phase space follows directly from the well-known Liouville's theorem for an 
N— particle system upon integrating over the six phase-space coordinates of N-l 
of these particles. 

The general approach may be formally understood as follows. The one- 
particle distribution function in kinetic theory has the physical meaning that 

f(v,x,t) dx dv (1) 


represents the (ensemble-average) number of atoms (or molecules) at time t in 

^ ^ I 

a spatial volume element dx about position x, with velocities dv about v. The 
time evolution of the one-particle distribution function is governed by the 
kinetic equation 
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where I 1 , which may in general depend on (v,x,t), is the force experienced by 
the gas atoms (or molecules) due to macroscopic internal and external sources 
and m is the mass of a gas atom (or molecule). Furthermore, (6f/6t) is the 
local rate of change of f(v,x,t) due to processes such as collisions that may 
add or remove the (ensemble-average) number of atoms (or molecules) from an 
element of phase space. For binary collisions in a gas, for example, one 
obtains the Boltzmann equation when the right-hand side of the kinetic 
equation is represented by the Boltzmann collision integral (a nonlinear and 
integral function of the distribution function). 

The time evolution or trajectory of a point in phase space (v,x) = 



determined by solving the equations of 
mass m 

_ dx(t ) 
vU) * d.' • 

?(v(t) ,x(t) , t) _ dv(t) 
m dt 

for this initial value problem. These 
the characteristic equations of the kinetic equation so that the left-hand 
side of the kinetic equation is the total time derivative of the distribution 
function. The kinetic equation thus reduces to 



There exists a class of physical problems to which the Io gas clouds may 
be related in which the gas density is sufficiently low that gas atoms (or 
molecules) along their phase space trajectories do not collide with each other 
but may collisionally interact with other material (such as ions and electrons 
in a rarefied plasma) or photons (i.e., photoionization or photodissociation 
processes) so as to change the initial ensemble-average number of atoms. In 
this case, the right-hand side of the kinetic equation is linear in the 
distribution function and may be expressed in the general form 

(ff) = “Vf + s > < 5 > 


motion for a gas atom (or molecule) of 

(3. a) 
(3.b) 

equations of motion are by construction 


where v and S are independent of f but may depend in general upon (v,x,t). 

Here v is an effective loss rate and S is an effective source rate for the 
ensemble-average number of atoms in an element of phase space. Using 
expression (5) in (4), the kinetic equation can be explicitly integrated along 
the phase space trajectory to determine the distribution function at time t: 


f (v(t),x(t),t) = f(v b ,x b ,t b ) P(t,t b ) 

+ f£ dt’ S(v(t’),x(t’),t’) P(t,t’) . 


( 6 ) 
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Here 


p (t 2 ,t 1 ) 



(7) 


is the probability that an ensemble-averaged number of gas atoms introduced at 
a time t^ exist at a later time t 2 because of the effective loss rate v. The 
distribution function at a time t > t^ is then given by the initial distribu- 
tion function reduced by this probability plus all new source contributions to 
the initial distribution function that have occurred in the time interval t^ 
to t with a proper decay loss included for these contributions. 

The description of the gas may therefore be constructed in terms of a 
large number of phase space trajectories along which the ensemble-averaged 
number of atoms evolve in time based upon a set of initial conditions for the 
distribution function which properly simulate the boundary conditions of a 
particular physical problem. More specifically, the description of the gas 
can be fundamentally reduced to describing the time evolution of an ensemble 
packet of gas along its phase space trajectory. An ensemble packet is defined 
as the infinitesimal volume element containing the ensemble-average number of 
atoms that initially began its trip along the trajectory in the time interval 
C b to bjj + dt. The physics of an ensemble packet appropriate for the lo 
neutral gas clouds is described in the next subsection. 

The basic principle of the general mathematical approach described above 
has been employed in planetary science to develop models for a number of 
studies involving collisionless or nearly collisionless atmospheres. The 
basic principle has been, in fact, used with increasing regularity in the last 
two decades. These studies have included cometary atmospheres, planetary 
atmospheres, and satellite atmospheres for which major references are given in 
Table 2. A more general review of exospheric theories is given by Fahr and 
Shizgal (1983). The studies for Io's neutral gas clouds are perhaps the most 
complex subset of these exospheric problems appropriate to the linear inter- 
action (5). This complexity results because of the inherent velocity, space. 


and time dependence of both the force F in the equations of motion and the 
effective loss rate v of neutrals in the plasma torus and also because of the 
form of the source term S which is described below. 

For the Io neutral gas cloud, the source term S in (5) is a spatially 
distributed source that describes the creation of secondary cloud atoms 


9 



produced by interactions of the raagnetospheric plasma with only the primary 
gas cloud atoms (or molecules) initially ejected from the satellite. For 
atomic gas clouds, the distributed source is caused primarily by charge 
exchange reactions between the extended cloud atoms and the corotating plasma 
torus ions. Elastic collisions between these ions and extended cloud atoms 
may also produce a contribution to the source term in (5) but this contribu- 
tion is usually significantly less important over the relative short lifetime 
of the cloud. Primary cloud atoms have velocities relative to Io that are 
generally slow compared to the corotational ion speed (~ 57 km sec *) and, 
when removed by the atom-ion charge exchange reactions, create fast secondary 
cloud atoms that have a rather tightly confined velocity dispersion centered 
on the corotational speed. The fast secondary cloud atoms are thus created in 
a velocity volume of phase space that is generally discontinuous with velocity 
volume of phase space containing the slow primary cloud atoms. For molecular 
gas clouds (e.g., SO2 or SO), electron-impact dissociative reactions very 
rapidly remove the primary cloud molecules (Shemansky 1984) and create other 
secondary cloud atoms (or molecules) that initially may have additional 
amounts of kinetic energy as dictated by the energetics of the reaction. 
Because of the rapid dissociative lifetime of molecules in the plasma torus, 
these dissociative sources are likely to be heavily confined near Io. A gas 
cloud of Io for the j neutral species may therefore be generally described by 
a primary and a secondary population of cloud atoms (or molecules) by formally 
dividing the one-particle distribution function (1) into two parts. 

f = f (primary) + f (secondary) , g v 

j j 3 

The kinetic equation (4) with the right hand side given by (5) may thus 
be formally subdivided for the j species into equations for the primary and 
secondary atom (or molecule) populations 


df (primary) 



_y f (primary) 
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d £ (secondary) 
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(9. a) 


(9.b) 


where the secondary source term may be written in the form 
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(primary) 


(9.c) 


where describes the creation rate of secondary cloud atoms (or molecules) 
of the j species from primary cloud atoms (or molecules) of the k species. 

The source term (9.c) therefore represents the secondary cloud atoms (or mole- 
cules) produced by all primary cloud atoms (or molecules). More generally, a 
source term similar to (9.c) may likewise be introduced if necessary to 
include third generation extended source terms produced by dissociation of 
secondary cloud molecules (e.g. SO S + 0). Limiting the present development 
to; second generation source terms, the description of a neutral cloud is then 
obtained by first solving (9. a) for the properties of the primary cloud atom 
(or molecule) population 


(primary) (primary) 

f j(v(t),x(t),t) = fj(v b ,x b ,t b ) Pj(t,t b ) 


(10. a) 


and then by secondly determining Sj and solving (9.b) for the properties of 
the secondary cloud atom (or molecule) population 


(secondary) 

f (v(t),x(t),t) = j dt" S (v(t"),x(t ,, ),t")P (t,t”) . (10. b) 

3 t 3 3 

s 

(secondary) 

where f t (v_ ,x„ ,t e ) = 0 by selecting the initial point of the secondary cloud 

j o 5 8 v 

phase space trajectory (v ,x ,t ) at a point on the primary cloud phase space 

s s s 

trajectory where the secondary cloud distribution originates and is hence zero 
by definition. 

In this paper, a general model for the primary cloud atoms (or molecules) 
is explicitly developed and will include expressions to calculate the relevant 
plasma Input rates supplied to the plasma torus by the primary cloud. A model 
for the secondary cloud atoms (or molecules) is also given but is more briefly 
outlined. For simplicity, the general model developed will adopt the language 
appropriate to a gas cloud of atoms but, with slight generalization, will 
equally well apply to a molecular gas cloud. For a gas cloud of atoms, the 
primary cloud population is much brighter than the secondary cloud population 



in the general vicinity of Io's orbit, but the secondary cloud population 
extends throughout and beyond the planetary magnetosphere and provides a 
background density and cloud brightness into which the primary cloud 
population merges. 

Physics of an Ensemble Packet 

A neutral gas cloud of Io may be described for the primary atom popula- 
tion discussed above by following the space— time history of neutral atoms 
ejected from the surface of the satellite (in low atmospheric density limit) 
or from the exosphere of the satellite (in the presence of a substantial 
atmosphere). Without loss of generality, a spherical exobase of radius R £ 

centered on Io is adopted. For a differential surface area element 
1 2 

dft ( = Rg sin$ d$dQ) of the exobase centered on the angular direction 
Q = (0,$), let the neutral flux (number of atoms ejected per unit area, per 
unit initial velocity interval, per unit time) in the initial velocity inter- 

> 4- a 

val a) to o> + da) and at the absolute time t^ be denoted by <j>(a>, ft, t^) . The rate 
at which these atoms are ejected from the surface element is then 

<|>(u), ft, t^) dft dm. (11. a) 

In terms of the one-particle distribution function of the preceeding subsec- 
tion, the flux <f> is defined by 

(primary) 

<K<d, ft, tjj) = m • ft f(u>, x^jt^) (11. b) 

where x^= (R^,, 0, *) is the spherical polar coordinate location of the surface 
area element dft. 

To describe the neutral cloud and its impact on the planetary magneto- 
sphere, it is sufficient to follow in time an ensemble packet of atoms for 
each neutral species initially ejected in an absolute time interval t^ to 
t^ + dt', where t' is defined at the flight time since the packet was initial- 
ly ejected from the exobase at time t^. For one such neutral species, the 
motion of a packet is depicted in Figure 1. A number of physical quantities 
will be defined for the packet in this section. The complete description of 
a neutral cloud for a particular physical quantity is obtained by integrating 
the packet contribution for this quantity over its complete flight time and 


over all initial exobase conditions. To develop this description, the 
location of the center of the packet at flight time t' may be denoted by 
rCt'j ItOjfi) in a coordinate system fixed ■ relative to the satellite. In the 
case of a sodium or potassium cloud, the additional force of solar radiation 
pressure experienced by the packet atoms as they resonantly scatter sunlight 
requires a more complex solution for the packet location (see Smyth 1979, 

1983) that depends in addition on the present value of the absolute time, 
t b + t’, and hence must be denoted by r(t',t b |m,n). This absolute time depen- 
dence is introduced because the magnitude of the solar radiation force is mod- 
ulated by the doppler motion that occurs between the sun and the packet atoms 
in the satellite ;frame as Io executes circular motion about Jupiter and as 
Jupiter moves on its elliptical orbit about the sun. 

As an ensemble packet of atoms moves through the circuraplanetary space , 
it may experience ionization, charge exchange, excitation, and even elastic 
collisions because of its interactions with the Io plasma torus. Neutral- 
neutral collisions are generally not important (except for recapture near Io 
and Jupiter where they are easily included) because the accumulated column 
density determined by the relatively short lifetime of cloud atoms is small. 
Electron impact ionization, photoionization, and charge exchange reactions 
remove atoms from the packet and hence determine its decay lifetime. Ioniza- 
tion processes also supply new ions to the planetary magnetosphere. Charge 
exchange processes modify the ion composition and energy of the Io plasma 
torus and also provide a flux of high-speed neutrals that (unless captured in 
the middle magnetosphere) escape from the planetary system. Electron impact 
excitation processes for neutral oxygen and sulfur atoms provide a photon 
emission mechanism for these two neutral gas clouds. Elastic ion-neutral 
collisions alter the trajectory of atoms in the packet and, if important, 
require that this momentum transfer information be properly incorporated in 
determining r(t’ , t b j to ) for the secondary atom population as described in the 
previous subsection. 

To describe the interactions of the ensemble packet of atoms with the 
plasma torus, it is necessary at an absolute time t to map a location R in the 
satellite fixed coordinate system where the packet location r(t',t b j«, ft) is 
specified to a location R in the magnetic (or, more correctly and generally, 
what shall be called the plasma) coordinate frame where the plasma properties 
are specified. The general transformation for this mapping, which is a 


function of the absolute time t, will be denoted as follows: 


R = T(t) • R + D(t) . (12) 

The tensor part of the transformation, T(t), includes the relative rotational 
motion and the tilt angle between the satellite coordinate frame and the 
plasma torus coordinate frame. The vector part of the transformation, $(t) , 
includes the more complex mapping required between the satellite frame and the 
plasma frame for a dipole offset of the origin of the magnetic and planetary 
spin coordinate frames as well as for the additional shift of the center of 
the plasma torus produced by an east-west electric field such as proposed by 
Barbosa and Kivelson (1983) and Ip and Goertz (1983). The location of a 
packet at flight time t' in the plasma coordinate system is then given as 
follows : 

r(t',t b |S,0) = T(t b + f) • *(t',t b |S,a) + fc(t b +t') . (13) 

+ A 

The velocity of the packet v(t ' , t^ |m ) at flight time t' in the plasma 
coordinate system is then determined directly by time differentiation of (13). 

In the remainder of this subsection, the ensemble packet contributions 

for a number of neutral cloud physical quantities of interest (e.g. density, 

photon emission rate, neutral loss rate, and various plasma input rates to the 

magnetosphere) are formulated. To facilitate this formulation, a hat " A " 

notation is adopted to indicate that a physical quantity is described in 

terms of its phase space representation (t',t. |u,£2) along its phase space 

b 

trajectory. Thus, for a general physical quantity Q, its ensemble packet 
contribution will be denoted by § = Q (t',t ju,S). For an observer in 
a coordinate frame (r^.t^), the space/time description of the physical 
quantity Q(rg,tg) is then determined, as discussed in a later section, 
from Q by evaluating its ensemble average <Q>. 

The number of atoms in an ensemble packet at flight time t' depend upon 
the loss processes suffered by the packet since its beginning at time t b « 

For an ensemble packet at a location r having a velocity v in the plasma 
coordinate frame, the lifetime of cloud atoms is denoted by x(r,v). The 
velocity dependence is required in the lifetime to calculate properly the 
cloud atom charge exchange reaction rates with the plasma torus ions. The 
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loss collision frequency for the packet is then given by 


v(t',t b |w,ft) = [T(r(t' ,t b |u,0),v(t' ,t b |w,fl))] 


-1 


( 14 ) 


and the probability that the ensemble packet of atoms at flight time t' exists 
is defined as follows: 


P(t',t b |5,n) = exp [-/|j v(t",t b |u,£2) dt"] . 


(15) 


The number of atoms in the packet at flight time t' is then given by 


n (t',t b |m,fi) d£2 dm dt' 


(16) 


where 


n (t',t b |u,fi) = P (t',t b |m,n)4(m,n,t b ) (17) 

is the number of atoms per unit surface element dfi at the initial surface 

location ft, per unit initial velocity volume dm at the initial velocity 10 , and 

per unit time interval dt' at the flight time t'. 

The number of photons emitted by an ensemble packet at flight time t' 

depends upon its photon excitation mechanism. The two important mechanisms, 

electron impact excitation in the plasma torus and solar resonance scattering, 

are described by an atom excitation rate (i.e. photons per atom per unit time) 

denoted by J,(r, v ) for an emission line with wavelength X. The position 
£ A s 

vector r of the packet in the plasma coordinate frame is required in 
calculating the electrom impact excitation rate for oxygen and sulfur atoms. 
The radial speed v g of the packet relative to the sun is required in 
calculating the solar resonance scattering excitation rate for sodium, 
potasium, oxygen and sulfur. The photon emission rate (i.e. number of photons 
emitted per unit time) at wavelength X for the packet at flight time t* is 
then given by 


(19) 


i x (t' ,t b |u,n) = J x (r(t' ,t b |w,fi), v 8 (f,t: b |J i n)) n(t' ,t b |u>,fi) 


The rate at which : neutral atoms are lost from the packet at flight time 
t' is determined from the time derivative of (17) and is given by 


9 (t 1 , tulw.fi) .a . 

( ) <|>(w,fi,t b ) dfi dw dt' 


( 20 ) 


which may be written in terms of a loss rate 


L(t' ,t b |u> ,fi ) dfi dm dt' 


( 21 ) 


where 

L(t',t b |u) ,fi) = v(t',t b |w,fi) n(t',t b Jw,fi) . (22) 

A 

Loss processes contributing to the collision frequency v in (22) occur for the 
first four types of neutral reactions summarized in Table 3, where the neutral 
cloud atoms are denoted by the general chemical symbol Y. (It should be noted 
in passing that to generalize the collision frequency v for molecular gas 
clouds, dissociative reactions would from this point forward have to be 
properly included.) It should be emphasized that the expression (21) is not , 
however, the net rate at which ions are deposited into the magnetosphere by 
the packet, since charge exchange reactions of type C in Table 3, while 
providing from the packet a new ion (Y ), simultaneously remove an old ion 
(X + ) from the magnetosphere by production of a fast escaping neutral. 

In general, to determine the impact of the neutral cloud on the 
magnetosphere it is necessary to simultaneously assess the contributions of 
the packets of several neutral clouds on the plasma torus because of reactions 
of type C in Table 3. It is of interest to determine these contributions not 

only for the net ion loading rate N but also for the plasma mass loading rate 

• • 

M, and for the ion energy input rate E. For a mixture of N neutral gas 
species ( j=l ,2, . . . ,N) , the instantaneous packet contributions to any one of 
these three input rates can be formulated (see Smyth and Shemansky 1983) and 
will be denoted by R(t',t b |m,fi) , where 
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R 


N 


li i L i 'Vi 1 


(23) 


Here Lj is the loss rate of the packet for the j-species defined in a parallel 
manner to (22), and is the production rate of the packet for the j-species 
as determined from the reactions of type C in Table 3. For the three cases of 

A. A. 

interest, and pj are defined in Table 4, where is the relative velocity 
between the neutral atom Y and the plasma coordinate frame at the time the 
neutral is lost from the packet, and is the thermal energy of the X + ion at 
the time it is converted to a fast neutral because of reactions of type C in 
Table 3. 

Combining the ionization reactions A and B in Table 3 into a single term, 
the collision frequency v j (t' .t^lai ,ft) for the j species may be divided into a 
contribution from ionization (I), charge exchange (CE), and charge transfer 
(CT) as follows: 


C(I) + c< CE > + 0< CT) 

j j j 


(24) 


The charge exchange collision frequency may be further divided into terras 


~ (CE) 

j 


N 


(t'.tjs.n) = i s5 k) (t- 


k=i J 



(25) 


A (k) 

that describe the individual charge exchange reactions. Here S. is the loss 
rate for the j-species neutral packet because of the charge exchange reaction 
that produces a fast k-species neutral (i.e. it is the product of the ion 
density for the ionized k-species times the reaction rate constant for the 
charge exchange reaction). Symmetric charge exchange (k=j) and nonsymmetric 
charge exchange (k+j) are included. The production rate P (t', t^Jra, ft) for 
all fast neutrals of the j species that are produced by all neutral packets is 
then given by 


P. 

J 




(26) 


so that the expression (23) may be rewritten as the sum of two different terras 


N „ „ 


R = •£, M V i 

J=1 J J 


(I) 


+ c <ci) r„. ♦ i I (si< k > 

1 ' J j*l k*l J J 1 J k 


A 

°k) 


(27) 
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The contribution of the charge exchange reactions to (27) is given by the 

^ A 

• • 

second term. This second term is zero for N, nonzero for M because of the 

A 

nonsymmetric reactions, and nonzero for E because of every charge exchange 
reaction. 

A 

In the case of N, the number of neutrals lost from the j-species packet, 

a 

A • 

L - , is divided between the net ion loading rate, N,, and the production of new 

a / or- 

fast neutrals, v j Wn j» escape from the system because of reactions of 

type C (Table 3): 


L. 

J 



+ v 


(CE)" 

. n . 
3 J 


(28) 


as expected on physical grounds. The fraction of L. that contributes to Nj 
is time dependent because of the oscillation of the plasma torus about the 
satellite plane and other more complex motions that are contained in the 
transformation (13). For M and E, the separate j-species contributions are 
coupled as described by (27). The overall mass and energy conservation is 
reflected in (23), where the first term is the contribution to the plasma 
torus by the neutrals lost from the packets, and the second term is the 
contribution produced by new neutrals being created by the packet-plasma 
reactions of type C in Table 3. 

A 

The cumulative value of the instantaneous rate R supplied to the plasma 
torus in a time period T by all of the neutral species packets, each of which 
has a current flight time t', may also equally well be formulated. For the 
absolute time interval from tg-T to tg, this cumulative value is given by 


[/^ R (t", t -t' |m,fi) dt"l dfi dm dt 1 
0 J 


where 


t 


£ 


( 


0 


t'-T 


if t' < T 
if t' > T 


(29) 


(30) 


The quantity (29) is useful in that, even though the instantaneous rate R 
integrated over all flight time is time dependent, its cumulative value (29) 
averaged overall flight time t' is not time dependent, the source and sink 
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processes for the neutral cloud are periodic with the period T. This will be 
discussed further in a later section. It should be noted here in passing, 
however, that in addition to the instantaneous contribution (27) or the 

^ A A 

cumulative contributions (29) of the neutral packets to N, M and E, the plasma 
torus may also contribute directly through reactions of type E, F and G in 
Table 3. These additional processes are not in the realm of the neutral cloud 
description and may be treated independently. 


Physics of an Ensemble Subpacket 

A neutral gas cloud of Io may be described for the secondary atom 
population by following the space-time history of all neutral atoms produced 
by the primary atom population. To accomplish this explicitly, the secondary 
atom population will be described in terms of an ensemble subpacket in an 
analogous way that the primary atom population is described in terms of an 
ensemble packet. The relationship between the packet and subpacket are 
illustrated in Figure 2. For a neutral gas subpacket of the j 1 "* 1 species, the 
number of atoms created with initial velocity g to g+dg along the primary 
packet trajectory r(t ' , t^ jm ,Q ) in the time interval t' to t'+dt' is given by 

Sj (t ' , t^ |td ,g ) dft dm dt' dg (31) 


where 


S. 

J 


(t',t b |«,o,|) = P (t 



(32) 


Here, g^ (g |m , t ' , t^ ) is the local initial velocity distribution for all 
secondary atoms created at time t' along this packet trajectory and is 
normalized to unity 


/ gj (!|m,a,f,t b ) dg = 1 (33) 

a 

and is given by (26) for the case of charge exchange creation of fast 
secondary neutrals. For a secondary population created by elastic ion-atom or 

A 

dissociative electron-molecule collisions, the expression for Pj must be 
appropriately generalized. The creation rate (32) in the subpacket formalism 
occupies a parallel position to the flux (ll.b) in the packet formalism. 


The location of the center of the subpacket, originally created at 
a flight time t"=0 with an initial velocity £ at the primary location 
r(t',t |w,fi) and corresponding to the primary trajectory time t^+t ' , will be 
denoted for t">0 by r(t",£; t ' , t b jw ,fi ) . In the plasma coordinate system, the 
location of the subpacket at a flight time t" is therefore given by 


r(t" ,t ;t ' ,t b |u> ) = ^Cth+t'+t”) • r(t",£; t 1 


* t b 


5,a) 


D(t b +t'+f 


(34) 


in an exactly parallel manner to the packet location (13). The velocity 
of the subpacket in the plasma coordinate system will be denoted by 
v(t",£; t',t b |w,fl) and is determined by the t" time derivative of (34). 

The ensemble subpacket contributions for different neutral cloud physical 
quantities of interest can be formulated as in the previous section for the 
packet ensemble. To facilitate this formulation, a double hat notation is 

adopted to indicate that a physical subpacket quantity is described in terms 

of its phase space representation (t“ t',t. |aj,n) along its phase space 

b th 

trajectory. The loss collision frequency for the subpacket of the j species 
is given by 


vj (t" ,| j f,t b |S,o) = [Tj(?(t" ,t; f,t b |S,o), ?(t",£; f,t b |S,n)) ] _1 (35) 


parallel to the packet expression (14) and the probability that the ensemble 
subpacket of atoms at flight time t" exist is defined by 

t" : + + *■ 

P^t',1; t',t |w,fi) = exp [-/ v (t"',5; t',t |u),0) dt”'] . (36) 

J f j b 

i. 

The number of atoms in the subpacket of the j 1 species at time t" is then 
given by 



t'.t^jw,!!) d SI dw dt ' d£ dt 


(37) 


where 


n_.(f,£; t’,t b |m,n) = (t • , t fa |u ,£2 ,£) Pj(t",£; t'.tjw.n) 


(38) 


In a similar manner to the packet quantities defined in the previous section, 
the subpacket photon emission .rate density l^(t",£; t',t b |u>,ft), the loss 
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rate density 1 (t",£; t',t^|u f fl), and the plasma input rate densities 
R (t",£; t',t^jio,fi) may be defined parallel to the expressions (18), (21) and 
(27), respectively. 


Kinetic Theory Distribution Function 

The kinetic theory distribution function for the neutral gas cloud of the 
jbk species as noted earlier may be divided into a contribution from the 
primary and secondary atom populations 

(primary) (secondary) 


_ .-f •f . -f v 

f.(v ,x ,t ) = f .(v ,x ,t ) + f.(v ,x ,t ) 
J o’ O O J o o o J o’ o* o 


(39) 


Here, v q and x q are the velocity and spatial coordinates of phase space for an 

observation coordinate frame, and t is the absolute time in the observation 

o 

frame. The primary and secondary contributions to (39) can be directly 
expressed, respectively, in terms of the primary packet density n and the 

A 

secondary subpacket density n defined in the previous two sections: 


(primary) 

VvvV = 


a + 

J dn /dm / n. (t’,tQ-t, |u,n) 6 (?q,£, ) 6(v 0 ,r, ) (40) 
exobase velocity flight J J 

space time 


(secondary) 

VvvV “ / * / ^ f dt ' J J dt " 

J exobase velocity flight velocity flight 

space time space time 

A 

In (40) and (41), $ and r are abbreviated notations for the ensemble 
packet spatial location and velocity 


nj(t",5; t’.tQ-t' |a),fi)6(r 0 ,rj)6(v o ,rj) 

(41) 


A 

r = r(t ' ,t 0 - t ' |a>,fl) , r = [^r r(t ' ,t b |u,fi)] t =( . _ fc , 
. b 0 


(42) 


and ^ and r are abbreviated notations for the ensemble subpacket location and 
velocity 


r = r(t",| ;t’,t 0 -t' |m,a) , r = [r(t",| ; t ’ ,t Q -t ’ |m,n) ] * (43) 


The quantities 6 ( r 0 .‘t j ) and 6(vQ,rj), respectively. 


cue uue 


3 J/O UAOA 


21 


projection operator and the velocity volume projection operator defined as 
follows: 


6 (r Q ,r) 


1 

A 3 r. 


1 if *o < *o < V' A? o 

0 otherwise 


5 (v Q ,r) 




1 if ^0 < ^0 < ^0 + A ^0 

0 otherwise 


(44) 


(45) 


Here R^ is the location of the ensemble packet (or subpacket) in the 
observation coordinate frame defined and mapped by an appropriately defined 
transformation of the form 




-► 

r 


+ D (t ) 
0 0 


(46) 


and = Rq is the velocity of the ensemble packet (or subpacket) in the 
observation coordinate frame. In addition, the quantities 4 3 r^ and A 3 v^ are, 
respectively, the spatial and velocity volume elements establishing the volume 
resolution element in the phase space of the observation coordinate system and 
correspond to the spatial volume Arg about Tq and the velocity volume Avq 


about v 


0 * 


Spatial Properties of the Neutral Cloud 

The spatial properties of the neutral cloud (e.g., the three-dimensional 
number density or the two-dimensional brightness of the cloud on the sky 
plane) as well as the physical ion input rates to the planetary magnetosphere 
can be determined directly from the kinetic theory distribution function (39) 
or alternatively from the expressions derived from the ensemble packet and 
subpacket in the previous sections. In the former case, one multiplies the 
distribution function by a physical quantity of interest, which may in general 
be a function of (^Q^^.t^), and integrates the product over the velocity 
coordinate space v^. In the latter case, one defines an ensemble average for 
the packet and subpacket physical quantity of interest, producing an identical 


result to the former method because of the nature of the spatial and velocity 
volume projection operators (44) and (45). The ensemble average approach will 
be explicitly developed below and is the description adopted for the numerical 
evaluation of the neutral cloud model discussed in the next section. 

a 

If an ensemble packet quantity (e.g., n, 1^ , N) and its subpacket 

A A A A ^ A 

quantity (e.g., n,I^,N) are denoted respectively as Q(t',t, |w,0) and 

A ^ A ^ 

Q(t" ,£; t',t^|m,£i) , and if the corresponding ensemble averaged physical 
quantity in the observer's coordinate frame at r^ and at an absolute time tg 
is denoted by Q(rg,tg), then 

A A 

Q(r 0 ,t 0 ) = <Q> + <Q> (47) 

I 

where 


A 

<Q> 


J da / dm 

exobase velocity 
space 


/” dt' Q(t , ,t () -t , Iu,n)6(r 0 ,r) 

flight 

time 


(48) 


A 

<Q> 


/ da 

exobase 


/ cEa 
velocity 


/ dt' 
0 


s D ace flight 

space tin)e 


/ d? / dt" 

velocity flight 
space time 


0(t",|; t',t 0 -t' |m,a)6(r 0 ,r) 


(49) 


where t and ^ are given by (42) and (43) and where the volume projection 
operator 6(rg,r) is defined by (44). To project an ensemble average physical 
quantity Q(Cq. t^) onto the observer's viewing plane, as is appropriate for 
the gas density to obtain a column density or for the volume photon emission 
rates to obtain an image of the cloud brightness, an integration along the 
line of sight in the observer's coordinate frame is required in addition to 
those specified in (48) and (49). For the simple case of no parallax with the 
Zg-axis along the line of sight, the column integrated quantity is given by 


Q(x 0’ y 0 t 0 ) = Q( ^0’ t 0 ) dz o* 


(50) 


In this manner, the basic properties of the cloud and also its interactions 
with the magnetosphere may be calculated. 



It should be noted that if the ensemble packet quantity and the packet 
trajectory are periodic with period T in their initial ejection time t^, 
that is, if 


A ^ A A A 

Q(t ' , t+T|u ,fi ) = Q(t',t|u,fi) 
r(t',t+T|u,fi) =r(t', tlu.fi) » 


(51a) 


(51b) 


it then follows that the corresponding ensemble average physical quantities of 
the neutral cloud (47), although time dependent in general, are also periodic 
with the period T 

Q(r 0 ,t 0 + T) = Q(r 0 .t 0 ) • (52) 

In general, the ensemble averaged value of the cumulative value [see (29) and 
(30)] of any packet quantity Q integrated over a time tg-T to tQ 


< [/£ Q(t",t 0 -t|u,fi) dt"] > 


is exactly equal to the ensemble average value of the instantaneous packet 

A 

quantity Q integrated over the time interval 


J <Q> dt* 

nn v 


as can be readily verified. A similar result also holds for the subpacket 
quantity. If (51) holds, the quantities (53) or (54) are also explicitly 
independent of the time tQ, as would, of course, be expected purely on 

A A 

physical grounds. For Q = L, this time-independence is clearly illustrated 
in the simple case of a constant flux 

4>(u,fi,t.)=<j)(u,fi) (51 


by the obvious result 


t n +T „ 
,->• r 0 ~ 


/ dr Q / t u <L> dt Q * = T / dfi / du <)>(u,fi) 

all 0 exobase velocity 

space space 


developed directly from (20) and (48). 
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The general periodic relationship (52) and the equality of (53) and (54) 
are useful in studying the impact of the neutral clouds on the magnetosphere. 

A period T of approximately 13 hours exists between Io and its location in the 
oscillating/rotating magnetosphere if east-west torus differences are assumed 
small or zero. If on the other hand, east-west torus differences cannot be 
ignored, a period T of approximately 32 days exists for recurrence of the same 
east-west position and System III magnetic longitude of Io. To illustrate the 
impact of the the neutral clouds, the packet expression (28) is summed over 
the various neutral cloud species and its ensemble average is taken to obtain 


^ N 

<L> =<&>+<£ 
3=1 


v . 
J 


(CE)~ . 

n . s 
J 


(57) 


This expression divides the instantaneous total packet neutral loss rate 
of the cloud <L> into its instantaneous net ion loading rate for the 


magnetosphere <N> and its instantaneous production rate of fast neutrals 
< Ij v j nj>. This partitioning is time as well as spatially dependent. 

The latter term of (57) is, of course, the rate at which subpacket fast 
neutrals are produced. If the packet periodicity (51) holds, the period 
averaged value of each term in (57) is independent of time, so that the 
average fractions of neutrals contributing to the net ion loading rate and to 
the production of fast neutrals are given, respectively, by 


*(N),+ x . v f A (lk ''(CT)^ s A A s 

f (r ) = < ) v; + v. n.> / < ) v.n.> 

° j=i J J 3 j=l J J 


and 


£ <fa St) t « 5®); , / < J c ; > 

o . , 3 j . , j j 
j=l J J j=l J 


(58a) 


(58b) 


where the brackets "<>" denote the packet ensemble average and the bar 
indicates periodic accumulation, that is 


t +T 
o 


t +T 
° nft 


<q> - <Q> dt; - / t Q(r o ,t;) dt; 


(58c) 


and where by construction 


f (K) ( x + f (£ast) ( . } , j 
o o 


(58d) 


4. Numerical Approach 

The ensemble average of a packet and subpacket quantity -given by 
(47) - (49) and its column integrated value given by (50) may be evaluated 
numerically. For the packet description, this is accomplished by dividing 
the exobase into many small area elements. On each area element, the 
initial velocity distribution as well as the absolute flux are specified 
independently. For each initial condition specified on an area element, a 
trajectory r(t' ,t b |w ,fi) for an ensemble packet is calculated for sufficiently 
long flight time t' to determine accurately the cloud properties. Along each 

A 

trajectory, the packet quantities Q of interest are calculated and accumulated 
in their appropriate spatial volume (or area of the sky plane) in such a way 
to determine the packet ensemble average (48) or its contribution to the 
column-integrated ensemble average (50). For the subpacket description, the 
source density rate S for the subpacket as defined by (31) is calculated in 
the packet description along the packet trajectory and has the analogous 
role that the flux (11a) has in the packet description. The source density 
rate S is numerically divided into many small packet flight time increments 
along the packet trajectory and from each increment the subpacket trajectories 
r(t",g; t',t b |(D,tt) are calculated for sufficiently long flight time t" to 
determine accurately the cloud properties. Along each trajectory, the 

A 

subpacket quantities Q of interest are calculated and accumulated in the 
appropriate spatial volume element of the observation coordinate frame. 

The packet trajectories r(t ' ,t b |m,ft) and subpacket trajectories 
r(t",f; t',t b |u>,f2) may be determined by solving the standard (Smyth and 
McElroy 1977) or modified (Smyth 1983) circular-restricted three-body 
equations of motion for each prescribed initial condition. The modified 
equations are required for the sodium and potassium cloud because of the 
additional effects of solar radiation pressure (Smyth 1979, 1983) on these 
atoms. The equations of motion are solved numerically by applying a straight- 
forward Runge Kutta technique. The time steps employed in the numerical 
solution must be selected small enough not only to determine an accurate 
trajectory but also to insure that the space and time resolution of 
the ensemble average quantities (47), which are quite sensitive to the 
oscillational/rotational motion of the plasma torus, are calculated to the 
level of desired accuracy. 


For a particular neutral cloud, calculation of the palcket quantity Q, 

A. 

the subpacket quantity Q and their ensemble averages (48) and (49 )- requires 
firstly that the relevant reactions and their rates (symbolically indicated by 
A-E in Table 3) and also the excitation rates (18) be identified and defined, 
and secondly that the space/time properties of the plasma torus be specified 
along the packet and subpacket trajectories. For the sodium cloud, this 
information is presented in the companion paper (Smyth and Combi 1986) in this 


issue. 


5 . Summary 

A brief review of Io’s local atmosphere and its extended gas clouds as 
well as their interactions with the planetary magnetosphere was given in 
section 2. The Io neutral clouds were seen to occupy a central position in 
understanding the coupled planet-satellite-magnetosphere system. Analysis of 
earth-based, rocket, or earth-orbiting satellite observations of neutral cloud 
emissions provides a direct, very promising, and still largely untapped 
approach for increasing our knowlege and understanding of complex physical 
phenomena in the system. Physical realistic models for the neutral clouds are 
essential in this analysis. These models must include the complex space-time 
interactions that occur between the neutral clouds and the plasma torus which 
have hitherto not been incorporated. The general Io neutral cloud model 
developed here includes these interactions and provides a general and solid 
base for future analysis of cloud data. 

A mathematical description for the general neutral cloud model was 
presented for Io in section 3. For a given neutral gas species, the 
description is divided into a primary cloud population that is initially 
ejected from the satellite exobase and a secondary cloud population that is 
created in a spatially extended volume by the primary cloud population through 
various time-dependent processes. For simplicity, the mathematical descrip- 
tion for the primary and secondary gas populations are explicitly developed 
for an atomic cloud, although they may be readily generalized for a molecular 
cloud. The description of the primary gas cloud may be fundamentally reduced 
to describing the time evolution of an ensemble packet of gas along its phase 
space trajectory. The description of the secondary gas cloud may likewise be 
reduced to describing the time evolution of an ensemble subpacket along its 
phase space trajectory. The physics of an ensemble packet and ensemble 
subpacket are presented and are also directly related to the kinetic theory 
distribution function for the gas cloud. The spatial properties of the 
neutral cloud as well as various ion input rates to the planetary magneto- 
sphere are obtained from their ensemble packet and subpacket values by 
evaluating their ensemble averages. A straightforward numerical approach for 
evaluating the model is outlined with appropriate references that define and 
present a method of solution for the equations of motion of the phase space 
trajectories for the ensemble packet and subpacket. 


As noted earlier, the general model is applied to the Io sodium cloud in 
a companion paper (Smyth and Combi 1986) in this issue. The application of 
the general model to other neutral gas clouds of Io (e.g., K, 0, S, S0 2 , etc.) 
will be considered in a series of future papers. 
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Table 2 


Some Studies for Exospheric Problems in Planetary Science 



Subject Reference 


1. Cometary Atmospheres Eddington 1910; Haser 1957, 1966; Keller 

1973; Keller and Thomas 1975; Keller and 
Meier 1976; Combi and Delsemme 1980; Combi 
1980; Festou 1981a,b; Combi and Smyth 1985; 
Combi, Stewart and Smyth 1986 


2. Planetary Atmospheres 

• Mercury Hartle 1971; Hartle, Ogilvie and Wu 1973; 

Hodges 1974; Hartle, Curtis and Thomas 1975; 
Smith et al. 1978; Curtis and Hartle 1978; 
Hodges 1980; Ip 1986; Smyth 1986 

• Venus Penner 1977; Hodges and Tinsley 1981, 1982; 

Rodriguez, Prather and McElroy 1984; Hodges 
and Tinsley 1986 

o Earth Chamberlain 1963; Hartle 1973; Bertaux and 

Blamont 1973; Fahr and Paul 1976; Penner 
1977; Hodges, Rohrbaugh and Tinsley 1981; 
Rodriguez 1983; Bishop 1985 


3. Satellite Atmospheres 

• Moon 

• Titan 

• Io 


Hodges 1973; Hartle and Thomas 1974; Hodges 
1980 

Fang, Smyth and McElroy 1976; Smyth and 
McElroy 1977; Smyth 1981; Shemansky, Smith, 
Smyth and Combi 1984; Hilton and Hunten 
1985; Ip 1985 

Carlson, Matson and Johnson 1975; Fang, 
Smyth and McElroy 1976; Smyth and McElroy 
1977, 1978; Matson et al. 1978; Carlson et 
al. 1978; Smyth 1979; Goldberg et al. 1980; 
Macy and Trafton 1980; Smyth 1983; Sieveka 
1983; Pilcher et al. 1984; Smyth and Combi 
1983, 1984, 1985; Sieveka and Johnson 1985; 
McGrath and Johnson 1986 
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Table 3 


Classif ication of Reactions Involving Neutral Species 


Neutral-Solar Reactions: 
A. Ionization 


Y + hv -*- Y + + e 


Neutral-Plasma Reactions: 

B. Ionization 

C. Charge Exchange 

D. Forward Charge Transfer 


Y + e ->■ Y + + 2e 

Y + X + ■> Y + + X 

Y + X** > Y* + X + 


Plasma-Plasma Reactions: 

E. Recombination 

F. Reverse Charge Transfer 

G. Dissociative Recombination 


Y* + e -*■ Y 
Y + + X + -► Y + X ++ 
XY + + e ->• X + Y 
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Figure 2. 


FIGURE CAPTIONS 

Motion of a Neutral Cloud Packet. The mathematical description 
for the primary population of the neutral cloud based upon the 
dynamic evolution of neutral cloud packets is depicted. The one 
neutral cloud packet shown on its trajectory at a flight time t'>0 
is ejected from a satellite exosphere area element dft with a 
velocity m at flight time t'=0. 

Motion of a Neutral Cloud Subpacket. The mathematical description 
of the secondary population of the neutral cloud based upon the 
dynamic evolution of neutral cloud subpackets is depicted. The 
one neutral cloud subpacket shown on its trajectory at flight time 
t">0 is ejected at flight time t"=0 with initial velocity f from 
the neutral packet. The location of the neutral packet at the 
later time (t">0) is also indicated. 
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